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Abstract 

66 ' 

We calculate the B B parameter, relevant for B —B° mixing, from a lat- 
I tice gauge theory simulation at (3 = 6.0. The bottom quarks are simulated 

|-y ^ | in the static theory, the light quarks with Wilson fermions. Improved smear- 

' ing functions produced by a variational technique, most, are used to reduce 

statistical errors and minimize excited-state contamination of the ground- 
state signal. We obtain l? B (4.33GeV) = 0.98^4 (statistical) l^g (systematic) 



> , 

\Q , which corresponds to B B = 1.4CCg (statistical)^^ (systematic) for the one- 

loop renormalization-scheme-independent parameter. The systematic errors 
' include the uncertainty due to alternative (less favored) treatments of the 

perturbatively-calculated mixing coefficients; this uncertainty is at least as 
. large as residual differences between Wilson-static and clover-static results. 

Our result agrees with extrapolations of results from relativistic (Wilson) 
heavy quark simulations. 
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I. INTRODUCTION 



The experimental observation of B —B° mixing allows, in principle, the extraction of 
the | Vtd | CKM matrix element The over-determination of the CKM matrix is a 

high-precision test of the standard model of particle physics and is regarded as a poten- 
tial harbinger of new physics. The dominant uncertainty in the extraction of | V t d | from 
experimental measurements is due to theoretical factors from non-perturbative QCD. The 
key factor is B B f B , where f B is the 5-meson semi-leptonic decay constant and B B is the 
"bag constant" for the 5-meson, defined as the ratio of the matrix element of the operator 
relevant for the mixing to its value in the vacuum-saturation approximation (VSA). 

There have been a large number of lattice gauge theory simulations which have calculated 
the f B decay constant; however, much less work has been done on the calculation of the 
B B parameter. The earliest result || suggested that the VSA works quite well; this result 
was unanticipated and is quite non-trivial, as was reiterated by Soni 0. Later results by 
other groups are surprisingly scattered, with significant disagreement in some cases || and 
with some results markedly different than that suggested by VSA. Here we argue that, in 
fact, most raw lattice data are consistent with VSA (including ours which are quite precise 
due to the use of improved smearing functions) and that groups differ due to their choices 
of how to relate these to the full-theory continuum value. We argue that although large 
systematic uncertainties remain due to unknown higher-order contributions in the mixing 
coefficients, it is possible to formulate the calculation in a way which is stable against changes 
in normalization (such as tadpole improvement). Our result is in accord with VSA and is 
also in agreement with the large-mass extrapolation of calculations || which use relativistic, 
rather than static, heavy quarks. 

Some of the first attempts at simulating the static theory calculated both the decay con- 
stant and the B B parameter [§,|J. However, the required perturbative matching coefficients 
were not known; these have since been computed by Flynn et al. [0. Their analysis showed 
that additional operators, not included in the first simulations, are required to estimate the 
B B parameter. 

Until recently, the simulation of the static theory was problematic because of excited- 
state contamination of the ground-state signal [jj^|10| . The development of variational 
techniques (|Tl],|r^] has finally allowed a reliable extraction of the decay constant. In this 
paper, we use a modern variational technique |Tj| to obtain accurate estimates of the lattice 
matrix elements and combine these with the mixing coefficients to calculate the static B B 



parameter. At two conferences [13,14], we have reported preliminary results for the value of 



B B from this simulation. 

Sec. [Tl] outlines the method of extracting the relevant matrix elements from lattice cor- 
relation functions; Sec. [Hj summarizes our numerical results. Sec. |IV| contains a summary 
of the perturbative-matching techniques which, rather explicitly, details our preferred way 
of organizing the calculation; we argue that our method reduces systematic errors in the 
matching coefficients which are then estimated in Sec. [VJ. In Sec. N\ a comparison is made 
to other groups as an illustration of the differences in the methods discussed in Sec. |V[ The 



conclusion follows as Sec. VII 



2 



II. NUMERICAL TECHNIQUES 



The static-light B B parameter is obtained from a combination of two- and three-point 
hadronic correlation functions. The required three-point function is 

C 3 , x (h, i 2 ) = EE(°| T f i) *(°> 0) X(*2, £ 2 )) I ) (2.1) 

which has a fermion operator inserted at the spacetime origin between two external B- 
meson operators, x- The times are restricted to the range t\ > > t 2 . We use the spatially 
extended 5-meson interpolating field 

X (x, t) = Yl /(r) f + f ) 7s &(*, x) (2.2) 

r 

where / is a smearing function chosen [12j to project out the ground state efficiently. The 
four-fermion operators, Ox (with X G {L, R, N, S}), are definedS as fTj: 

0l = &t„(i - 75)^(1 - 75)? 

Cfl = 67^(1 + 75)g&7 M (l + 7 5 )g 

0* = (26(1 - 7 5 )g6(l + 75 )<? + 26(1 + 75 )g6(l - 75 )g 

+67^(1 - 75)967^(1 + 7 5 )g + 67^(1 + 7s)gH t (l - 7 5 )g) 
Os = 6(l-7fe)g6(l- T5 )g (2.3) 

The operators Or and are introduced in the lattice and contribute towards Ol because of 
the poor chiral behavior of Wilson quarks. The operator Os is introduced in the continuum 
and contributes because of the matching of full QCD to the static theory. 
With the smeared-sink-local-source (SL) two-point function defined as 

C 2 (ti) = E(° I 7Xx(*i,£i) K0,0) 7475 9(0,0)) I 0) (2.4) 

the "raw" lattice-static parameters, B x , are calculated via the ratio of three- and two-point 
functions: 

Bx{tlM) ~ lc 2 ( tl) c 2 (t 2 ) > (2 - 5) 

The B B parameter itself can then be determined from the Bx = Bq x , extracted from fits 
of the Monte Carlo data to the form of Eq. (|2.5[), as the linear combination 



B B = Z Bl B l + Z Br B r + Z Bn B n + Z Bs B s (2.6) 



1 We choose a standard normalization for which the VSA value for Ol is (8/3) f B m B . 
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where the perturbatively-calculated mixing coefficients, Zb x , are defined in Sec. [IV[ Rather 
than this "fit-then-combine" method, our quoted results will be from the "combine-then-fit" 
method: 

B B (t 1 ,t 2 )= Yl Z Bx B x (t 1 ,t 2 ) B B (2.7) 

X=L,R,N,S 

For infinite statistics, the two methods should give identical results. 

We exploit time-reversal symmetry by averaging the correlators over t and T — t, where 
T is the length of the lattice in the time direction. We fix one of the times, t%, in Eqs. ( ^.5| ) 
and ( |2.7| ) and vary the other, t 2 \ the result is fitted to a constant. The fits include correlations 
in time, but not in the chiral extrapolation (a choice forced upon us by our limited statistics). 
The entire fitting procedure is bootstrapped (see, for example, Ref. [nj) to provide robust 
estimates of the statistical errors. An estimate of the systematic error due to the choice of 
interval is made by calculating the variance of the results from using all "reasonable" time 
intervals around our favorite one. 

A major problem with simulations that include static quarks is that the signal-to-noise 



ratio decreases very quickly with time |]8|, p!6|JT7[| ; therefore, the operator which creates the 



5-meson must project onto the ground state at very early times — before the signal is lost 
in the noise. Experience with the calculation of the f B decay constant in the static theory 
has shown that reliable results can be obtained only if the 5-meson operator is smeared with 
a very accurate "wave function," which can be obtained from a variational calculation on 
the lattice. We use the same smearing function as was used in our calculation of f B in the 



static approximation. This was obtained from the variational technique, called most ||T2 
which we have developed for this purpose. 

To demonstrate the effectiveness of the smearing function produced by most, we show 
in Fig. [I] the effective-mass plot (In C 2 S (t) / C 2 (t + 1) versus t + 1/2) for the two-point 
correlation function using a local (delta-function) sink at time t and an optimally-smeared 
source at time 0. The effective-mass plot has plateaued at small t (indicating the absence 
of significant excited-state contamination) before the signal-to-noise ratio has degenerated, 
so that a very precise mass and amplitude can be obtained by fitting over an early time 
range. If, instead, the same smearing function is used at the sink, with a local (delta- 
function) source, then it will still effectively remove excited-state contamination. Yet, as 
demonstrated in Fig. |^, this fact is obscured by much larger statistical fluctuations. (Since 
the spatial points are summed over at the sink to project out zero momentum regardless of 
which smearing function is used, smearing at the sink provides only marginal improvement 
in the signal and increases noise. In contrast, smearing at the source greatly enhances the 
signal and decreases the noise. For the local source the static quark is restricted to the 
spatial origin, and thus the statistics are poorer [f|.) 

We note that once an "optimal" smearing source has been obtained from the two-point 
function using a variational technique, it can be used directly in other calculations. The 
three-point function does not need to be formulated as a variational problem, although 
ground-state dominance should still be monitored using the mass splitting between the 
excited and ground states. 

The static quark never evolves in space from the origin because the four-fermion operator 
is at the spacetime origin. The 5-meson operator is constructed by smearing the light 
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quark relative to the heavy quark (Eq. (|2.4j) ). Fig. ^ shows a schematic of the quark flow 
resulting from the Wick contraction of Eq. (|2.1| ). The resulting two-point correlators are 
smeared-sink-local-source (SL) correlators, which are much more noisy than local-smeared 
(LS) correlators (as argued above) even though in the infinite-statistics limit they are equal. 
Since three-point functions are, in general, noisier than two-point functions, the "effective- 
mass" plots for these are even noisier than that for the SL two-point function; it would 
be hopeless to get a precise result for a static-light-meson matrix element without using a 
prohibitively large number of configurations. But fortunately, because the B B parameter is a 
ratio of matrix elements (Eq. ( |2.5| )), the noise is reduced due to the cancelation of correlated 
fluctuations between the numerator and denominator. 

It has been argued that the product B B f B and perturbative corrections to it should be 
calculated directly since it, rather than B B , is the phenomeno logically- important quantity. 
But there are several compelling reasons for calculating f B and B B separately. Firstly, 
although the calculation of B B , as for B B f B , is intrinsically more involved than is that of f B 
(both analytically, in the determination of perturbative corrections, and computationally), 
the numerical value of B B is more stable than is the value of either B B f B or f B . Certainly, 
f B is a very important physical quantity in its own right; it should be and is calculated 
separately. For this, one need only calculate a two-point function. However, the statistical 
fluctuations for f B are quite large; without the use of a reliable smearing function obtained 
variationally, excited-state contamination can be substantial and can mislead interpretation. 
(This may explain the scatter in the world summary of lattice calculations of f B |9|,|10|j.) 
Also, since its lattice-spacing dependence is rather large, especially when using the static 
approximation, its continuum extrapolation is delicate and prone to large systematic errors. 
Much computing effort is required to evaluate this simple quantity. However, B B (or B B f B ) 
requires the calculation of a three-point, in addition to a two-point, correlation function. 
Since it is more involved, it is usually determined in a secondary calculation after the primary 
calculation of f B and so fewer groups are likely to calculate it. Yet, as borne out by our 
data, since B B can be extracted from a ratio of three- to two-point functions which are 
strongly correlated, a quite precise value can be obtained, with an optimal choice of smearing 
function, from relatively few configurations. The calculational overhead (both computational 
and analytical) is large compared to the computational expense. Thus a handful of groups 
can fix precisely the value of B B once and for all, leaving for the wider community the task of 
applying improvements in algorithms and computers to the simpler f B . In the future, B B (in 
contrast to f B and B B f B ) need not be recalculated with every generation of improvements. 

Secondly, just as the numerical value of B B is stable because of cancelations of corre- 
lated fluctuations in numerator and denominator, we argue that so too are its perturbative 
corrections when linearized as is demonstrated in Sections |V| and |V]. The perturbatively- 
calculated coefficients for B B are likely more reliable than those for the product B B f B . 
Likewise, these are less likely to need updating with the next generation of improvements in 
analytic methods. 

Thirdly, it seems as though VSA is a surprisingly good approximation for the B B pa- 
rameter. This is an important qualitative statement, of use to model builders, which should 
not be obscured by poor-statistics attempts to calculate the product B B f B . 
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III. NUMERICAL RESULTS 



The simulations were carried out on a 20 3 x 30 lattice, calculated on 32 gauge config- 
urations, at (3 — 6.0. (This number of configurations is more than adequate for a precise 
estimate of the B B parameter with small statistical error since an efficient smearing function 
is used. The use of an ad hoc smearing function would have required an order of magnitude 
more configurations.) The simulations were quenched; the gauge configurations were gener- 
ated using the standard Wilson pure-glue action. The gauge configurations were fixed into 
Coulomb gauge. (An ultra-conservative gauge-fixing convergence criterion was used such 
that V-A was decreased to less than 10~ 9 its unfixed value.) The gauge-fixing was done 
only to choose smearing functions, but since these cancel in ratios of correlation functions 
all results are gauge-invariant (in the infinite-statistics limit). Wilson light-quark propaga- 
tors, with hopping-parameter values, k = 0.152, 0.154, 0.155, and 0.156, were used in our 
analysis. The value of kappa-critical used was 0.157 and the value of kappa-strange was 
0.155 @. 

Fig. |] shows, for the operator Ol, the ratio of the three- and two-point correlation 
functions .Bl^i,^), Eq. (|2.5| ), which asymptotically equals Bl for large Euclidean times. 
(In the figure Bi(ti, t 2 ) is graphed as a function of t 2 with t\ = —2 held fixed.) In fact, "large 
times" are remarkably small (>~ 2) because of the effectiveness of the smearing function in 
efficiently eliminating excited-state contamination, a fact supported by Fig. [I]. 

As with any lattice calculation of correlation functions, there is freedom in the choice of 
fit range and a balance needs to be struck between fitting over too-early times, for which 
systematic errors due to excited-state contamination may be non-negligible, and over too- 
late times, for which statistical errors will be unnecessarily large. In Fig. [| we display a 
t min -plot: the values for the fits of the raw B^ value (at k = 0.156) plotted for several 
choices of fit range. (All of our fits take into account the correlations in Euclidean time 
using the full-covariance matrix. For our central fit range, the values of the fits differ little 
whether or not the correlations are included.) The flatness of the plateau in Fig. |] reflects 
the insensitivity of the fitted value to the choice of fit range. For this and other plots 
we choose as our central values ti = —2 and 3 < t 2 < 6, a moderately-aggressive choice 
which has good x 2 /dof (0.83/3, 0.59/3, 0.41/3, 0.68/3 for k = 0.152, 0.154, 0.155, 0.156, 
respectively), small statistical errors, and fit-range systematic errors which are smaller than, 
but comparable to, the statistical errors. The fit-range systematic errors are determined 
from the standard deviation of all "reasonable choices." 

Figs. H an d | show similar plots for the raw lattice values for B R , B N and B$, 
respectively. Fig. |^ shows the ratio of correlation functions defined in Eq. ( ^.7[ ) from which 
the desired B B parameter is extracted. Again, the plot plateaus early with small statistical 
errors. Fig. pi] shows that, again, the value is insensitive to the choice of fit region. (For 
our central choice of fit range, the x 2 /dof are 0.74/3, 0.57/3, 0.43/3, 0.67/3 for k = 0.152, 
0.154, 0.155, 0.156, respectively.) We could also calculate the final B B parameter from 
the appropriate linear combination of the four fitted raw values B^, B B , B^ and B$, as in 
Eq. (pD. The x 2 /dof are good for B R (0.71, 0.67, 0.55, 0.33) and B N (0.60, 0.42, 0.36, 0.40). 
The worst x 2 /dof are for B$ (1.38, 1.30, 1.14, 0.85). The two procedures, combine-then-fit 
versus fit-then-combine, could give different answers in principle (for finite statistics), but 
in practice we see little difference. 
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K = U.lOZ 


K = U.154 


At = (J. 155 


k = U.lob 


K c = 0.157 


B L 
Br 
Bn 
-IBs 


1.01^(1) 
0.96iJ(l) 
0.97±1(3) 
1.00^(2) 


1.021^(1) 
0.96^(2) 
0.96l| (4) 
1-0011(2) 


1.02±|(1) 
0.95±l(2) 
0.96±|(4) 
l-00±l(3) 


1.03±g(2) 
0.95±§(2) 
0.96li(4) 
1.01±|(3) 


l-03±g(2) 
0-95lg(2) 
0.95l|(5) 
1.0ll|(3) 


B B (m* b ) 


0.95^(1) 


0.961^(2) 


0.961^(2) 


0.98l|(2) 


0.981^(3) 



TABLE I. The raw lattice B parameters for the operators Ol, Or, On, and Os which appear 
in the lattice-continuum matching, and the linear combination Bq{uii = B B as a function of k, 
and extrapolated to k c . The first errors are statistical (bootstrap) and the second are systematic 
due to choice of fit range. Note that Os has a VSA value different from that of Ol', with our 
normalization for the raw B parameters (a common denominator equal to the VSA value of Ol) 
—^Bs would identically equal 1 if VSA were exact. 

As shown in Table [L each raw lattice B parameter is close to 1.0 with small statistical 
errors, so our final valuep for B B {m\) is also close to 1.0, the VSA value, with similarly small 
statistical errors. 



IV. PERTURB ATIVE MATCHING 



To calculate the continuum value of the B B parameter, our "raw" lattice results, listed 
in Table |, must be multiplied by a lattice-to-continuum perturbative matching coefficient. 



After we finished the first analysis of our data [13[], we found that our value for B B was 
approximately 30% higher than the result of a similar simulation by the UKQCD collabora- 
tion ||21|| . We suspected that this difference was due to more than just the difference in the 



actions. This motivated us to do a very careful study of the perturbative matching, using 
the results in the literature, so that we obtained the "best value" of B B using the informa- 
tion available to us. (This is discussed further in Sees. 0, [VI], and jVIIj .) We also studied 



the systematic errors in the perturbative matching to find the reason for the disagreement 
between UKQCD's result and ours. 

For convenience, we shall refer to the AB = 2 effective Hamiltonian, obtained from the 
standard model by integrating out the top quark and the heavy vector gauge bosons, as the 
"full" theory although this is also an effective field theory. The perturbative matching is 
broken into two stages: full QCD to the continuum-static theory and the continuum-static 
theory to the lattice-static theory. For the matching of full QCD to the continuum-static 
theory, the relevant perturbative results have been calculated to do a next-to-leading-order 
analysis of the log(/i/mb) terms. The use of renormalization-group-improved perturbation 



2 



Bb is evaluated at m£, which is the scale at which the running mass is m{m^) 



in 



4.72 GeV 0JT|,|O| 
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theory reduces the renormalization-scheme dependence and the effects of the different ways of 
defining 75 in dimensional regularization p2| . Two scales are necessary for the perturbative 
matching: the scale, \i b = 0(m 6 ), of the matching to the full theory (we choose = m£, 
where ml is defined as mentioned earlier in footnote an d the scale, fi, of the matching 
to the lattice theory (we choose /1 = q*, which is determined from the Lepage-Mackenzie 
scale formulation p3] as discussed later in this section). Also, as emphasized by Ciuchini et 
al. []24| ], it is important to check the stability of the perturbative coefficient at next-to- leading 
order as the renormalization scale is changed. 

We choose, as do others pT] , pp] , pB| , to evaluate the full-theory operator, , at \i h : 

i=L,S 

where terms of order 1/m have been dropped. We use a double-argument notation similar 
to Ref. [26] to emphasize that this matching of the continuum-static theory to the full theory 
involves two theories (/ and c) and two scales ([i b and /i). C^(fi b ; /i) includes a running of 
the scale in the continuum-static theory which can be written explicitly due to the form 
of the solution to the renormalization group equation (RGE) for the coefficients (see, for 
example, pTfl ). 

C f ifa\ $ = C f j(fi b ; {r g expj - J ^ 

= C% b ;fi b )(u T y..( f i b ,fi) (4.2) 

Since we focus on the transformation of the operators, we treat the coefficients, C, as 
a row vector and transpose (T) the matrix U to be consistent with the common no- 
tation for U p8j|29| which treats the coefficients as a column vector: [C T ) c /{^\ = 
U?fji, fi b ) {C T )f{fi b ; fi b ) for which 

The superscript-label c indicates that the variables are for the continuum-static theory in 
which some degrees of freedom have been removed. Notice that the continuum-static scale- 
evolution matrix scales only the static-theory argument of the coefficient. Thus, Eq. ( |4.1|) 
becomes 

(0 / W>= E C% b -^ b ) (C/ T )^ 6 ,/i) (O^)) (4.4) 

i,j=L,S 

which is read, right-to- left, as "The static theory operator is scaled from \i to \i b where it is 
matched to the full theory." 

An alternative, not used here, is to evaluate the full theory operator at the same scale \x 
as is the static-theory operator, so that 

o'W) = E (om 

i=L,S 

= E (u T )%,n b ) C f ^ b]f i b ) (u T ) C ^ b ,fi) (4.5) 

i,j=L,S 
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(The generalization to multiple full-theory operators would include full-theory subscripts on 
0?, {U T y ', and C^.) Eq. (|4.5| ) reads, right-to-left, "The continuum-static theory operator 
is scaled in the static theory from \x to ^ where it is matched the full theory and then scaled 
back from //& to \i in the full theory." If U is treated to lowest order, summing neither the 
leading nor sub-leading order logarithms, then this reduces to the approach used by Flynn 
et al. who do not use the RG. The full-theory anomalous dimension appears there since 
this approach includes running the scale in the full theory. 

Returning to Eqs. ( |4.1| ) and (|4.4| ), matching in the continuum (with ^ = m b and H = q*) 
gives 

O f (mt)) = C?« q*) <Pi(q*)) + K; q*) (O c s (q*)) (4.6) 

We use the solution of the RG equation for a matrix of operators which is discussed by 
Ciuchini et al. and Buchalla in more detail. 



CfK;g*) =C f L c (mt;mt) 



a c s {ml)\ Po - L \ a c s (mt) - a c s (g* 




47T 



-PIl 



C f s c (mt; q*) 



+C ; s c (mt;m 



a c s (g* 



<x c a (q* 



with 



Po,i = (iQiJ (26q)) and p lti = [p 0>i (lujlo^ - &1/&0) 



(4.7a) 
(4.7b) 

(4.8) 



In Table [TI| we list the values of the anomalous dimensions of the various operators required 
in this calculation (all calculated using the naive dimensional regularization scheme). The 
coefficients from the first and second terms of the /5-function are defined as 
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Air 



47T 



16tt 2 



102 ~ fnj 
16vr 2 



(4.9) 



To obtain the leading-log expressions from the explicit solutions of the renormalization 
group equations that we quote, all quantities with a subscript 1 are omitted. In Eq. (|4.7b[) 
the higher-order terms of U c have been dropped when multiplied by C s c because Cg° is 
of order a s . We found that the inclusion of the Cg{U T )° s L term in our analysis was less 
than 0.05% of the C[ C (U T ) C L L term; this is smaller than the few percent effect which was 
quoted in Refs. p5"lj26l . Our ratio of the coupling at /i to that at fib was close to 1 because the 
automatic scale-setting procedure selected a scale q* which was close to m£. As the difference 
between the scales \i and ^ gets bigger, (U T ) C S L , which includes the leading off-diagonal 
terms in the anomalous dimension matrix, gets larger. 

We will now discuss the matching of the continuum-static theory to the lattice-static 
theory. The relevant perturbative calculations have been done by Flynn et al. 0. We want 
to calculate the full theory at m£: 



9 



Ref. 



their 
notation 



our 
notation 



value 



in 
(i) 



Ciuchini et al 
Gimenez [30 



1 23], Buchalla pi 



47- 



MS 



Ciuchini et al. [25], Buchalla [26 



Gimenez [30 



7 17 
7ll 

-Id ( 2 ) 

16 j\' 



7i, 



721 
(0) 
722 

do 
-J 



MS 



202 + 26^ 



16nj 



Ciuchini et al. 
Ciuchini et al. 

Ciuchini et al. 

Ciuchini et al. 

Ciuchini et al. 



|2J, Buchalla p6| 
|2§, Buchalla |26| 

||H, Buchalla |6| 

|| ; Buchalla g§ 

|25|, Buchalla p| 



Po,s 

Pl,L 



4 
3 
_8 
3 

26 



To 



s,s 



P0,L 



2b 



L.L 



L.L 



Duncan et al. [11 ,31 , 32, 33, 3 



7o 
7i 



lo.A 



' 1l,l 

P0,A 
f 

PlA 



-251 
9 



-4 

56tt 2 
27 



+ 



20n/ 



Buras et al. 35 



Buras et al. [35] 



-7 + f»/ 





TABLE II. Anomalous dimensions as defined by various groups and used here. The p's are 
defined in Eq. (|4.8| ) . All the results have been calculated using the naive dimensional regularization 
scheme. 
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C 



fc 



m 



s y m b^ 



O l L (a) 
O l s (a) 



+ Z 



O l R (a) 



+ 



O l N (a), 



+ 

(4.10) 



where Z (g*; a) relates the bare lattice-static theory matrix element to the renormalized 
continuum-static theory matrix element. After linearizing the product C^ c {m\\ q*) Z cl (q*; a) 
and allowing a separate coupling for continuum-static (a£) and for lattice-static 



a s ) we find: 



m. 



«sK) \ °- L L + «sK) - <* c sk* 

v a c s (q*) J { An 



a l s (q A 



(41n(g*V) - 21.7) J 



1 ( fa c s {mt)Y°> L 



a c s (q* 



An 

( a c s (mi) Y°' s 



An 



O l L (a] 



+ 



+ 







r 


















4tt 



-1.61) (^(a)) + ^(-14.4) (^(a) 



8) tefa^ 



(4.11) 



where we have updated the results of Flynn et al. @ by including {U T ) C S L |25|j26|] , by choosing 
the convention that the static-light two-point function be fit to the Ae~ mt model 



21.36 



and by including tadpole improvement [37 



Throughout this paper we assume the convention that the f B decay constant is extracted 
from the heavy-light correlators using the model Ae~ mt . Using this model changes the heavy- 
quark wave-function renormalization integral, denoted e, to a reduced value (see Eichten 
and Hill |36[). As mentioned by the UKQCD collaboration [21], this changes the Dl=— 65.5 
of Flynn et al. (the additive constant in the matching of the continuum-static O c L operator 



to the lattice operator) to D 



(R)_ 



-38.9. However, also appears in Z c \\ thus, the final 



values for the coefficients of the B parameters are independent of this choice if the ratio 
is linearized in a c and a 1 . In addition, any tadpole-improvement effects alter the three- 
point function by twice as much as each two-point function; linearizing the ratio cancels 
these effects exactly. However, when considering the three-point function and two-point 
function separately, one ought to include the effects of tadpole improvement. This changes 
the — — 38.9 to Dj^' tad ——21.7, as in Eq. ( f4.11|) . The large perturbative factors of the 



wave-function renormalization largely cancel in the expression for the B B parameters. 

To calculate the coefficients of B 
the static approximation is required 



To calculate the coefficients of B B} the renormalization coefficient of the axial current in 



52],|5B|,|54],|56],p8| ; we linearized the results quoted in 



Duncan et al. [ p 



Z cl 



A 



An 



-pU + 



An 



An 



(4.12) 

(21n(<fV) - 18.59)) 
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where the —18.59 is from using the mentioned above as well as including tadpole im- 
provement. If tadpole improvement had not been used, then this value would be —27.16. If 
e had been used instead of e^, then this value would be —40.44. As long as one is consistent 
between Eqs. ( 4.11 ) and ( |4.12| ), these effects cancel out of the linearized result for B f 



'B- 



The perturbative coefficients for the B B parameter can be obtained by dividing the 
four-fermion results by the square of Za and expanding the expressions linearly in a s . 



B 



B 



, <x c a {<? 



P0,L- 2 P0,A 



1 + a »<\-"-^ ( pU - 2 PU ) + 



47T 

+ ^1 (15.41) 



<x B \ m b 

47T 



26 

y 



47T 



+ 



+ 



a c s (q*) J 
' a c s {mt)\ pC °< L ~ 2pC °' A 



a c s (q* 



Atx 



B, 



' a c s (mt) \ pC °< s ~ 2pC °< A a c s {mt) 



-lM)B R+ ^(-UA)B N 



Itt 



47T 



B 



B 



Zb l Bl + Zb r Br 



Zn„B 



N 



-8)B, 



Zb s Bs 



(4.13a) 



(4.13b) 



where Z Bx = Un{Z x /Z 2 A ) = Lin ((C% Z$) / (C{ c Z%) 2 ) , X is one of {L, R, N, S}, and "Lin" 
signifies that the ratio is linearized as explained later in Sec. 0. The wave-function normal- 
ization factors of the quarks cancel between the numerator and denominator; no "tadpole" 
factors are required for this calculation if the coefficients are linearized. We also note from 
the values in Table |TJ that Pq L — %Pq a * s identically zero. However, Pq S — 2p^ A and p\ L — 2p\ A 
are not. If Eq. ( |4.13b| ) were expanded into explicit ln(///mj,) terms, then to first order the 
perturbative matching coefficients would not contain any logs. 

To calculate numerical values of the coefficients, we choose values for the scales /x& and fi, 
and for the couplings a c s and a l s . For a l s , we use a v , the coupling introduced by Lepage and 
Mackenzie P3"| . We use the plaquette value — lnWn = 0.5214 at (3 = 6.0. In the quenched 
approximation (n/ = 0), 



- \n(W u 



Air 



-a 



v 



3.41 



I — ax 



3.41 



(1.19) 



(4.14) 



which uses a lattice coupling which evolves with the form 



(4.15) 



where the j3 are defined in Eq. ( |L9|) . Eq. (|4.14|) defines a v and gives A v a = 0.169. 

Because the continuum-to-lattice matching is known only to one loop, these perturbative 
expressions are sensitive in principle to the value of the scale used in the matching. This 
dependence can only be reduced by calculating higher-order loops. However, Lepage and 
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Mackenzie |23| have described a plausible procedure for determining the scale and they have 



successfully tested this method for a number of quantities. 
The Lepage-Mackenzie scale q* is obtained from 

2 \ Jd 4 qf(q) \n(q 2 ) 

lnM > = WM (416) 

g*o = expQ (ln( ga ) 2 )) (4.17) 

where /(g) is the finite integrand of the lattice graphs; note that f(q) is defined by assuming 
that all the perturbative expressions are expanded linearly in the coupling. We used the 
integrands of Flynn et al. 0. (These have been confirmed by Borrelli and Pittori f39| .) In 



Table 01, we show the value of the scale for several operators. Our value for the scale for the 



static-light axial current, q*a = 2.18, agrees with the calculation by Hernandez and Hill |37 

The Lepage-Mackenzie scales for the individual operators in Table [H"T| are all around 2.0; 
however, the combined operator for B iuii has a lower scale of 1.22. (The scale quoted for 



B taa in the original preprint and conference proceeding was incorrect.) Morningstar [4(| 



also found very low scales for the perturbative renormalization of the quark mass in NRQCD 
(also see the comments by Sloan ^TJ); though this could be related to renormalon effects j42| . 



Using the scale of 1.22 gave large perturbative corrections. The Lepage-Mackenzie scale- 
setting procedure could be confused by taking the ratio of matrix elements of two operators 
that are approximately the same (obviously it would be inappropriate for the case of two 
equal operators because / would be identically zero). We chose to use the scale of 2.18/a as 
this is a typical scale for both and 0{^ 11 . 

We used A§£ D =0.175 GeV from D uncan et al. fLl|. They chose values for a 1 obtained 



from the charmonium system due to the low systematic errors. Although they do not quote 
a value for a~ l at /?=6.0, they did extrapolate A v from a -1 at (3=5.7, 5.9, and 6.1 in order 
to find a~ l at /3=6.3. We used this idea to interpolate to a~ 1 =2.1 GeV for /3=6.0. We also 
used their method for calculating m£; however, our number differs slightly because of the 
difference between the form of Eq. ( [4.15| ) and 



P 2 o ln(£ 



(4.18) 



With the full-to-continuum scale set as //& = m£ = 4.33 GeV and the continuum-to-lattice 
scale set by \ia = q*a = 2.18, we find a^(m£) = 0.21 and a l s {q*) = 0.18. 

Using a Monte Carlo technique, we estimated the error on the static B B parameter due to 
varying the values of the parameters used in the perturbation theory. A sample of one thou- 
sand was generated using uniform deviates for the renormalization scale, lattice spacing, the 
continuum Aq CD , and the bottom quark mass. The central value for each "input" -parameter 
distribution was set equal to our best value, based on those used in references |p|^7| . Rather 



than assume that the input parameters are known to three significant figures, we took up to 
20% of this value to be the standard deviation for each input parameter. The final results 
were sorted numerically and the 68% error range was taken as the "output" error. This 
procedure should produce more accurate estimates of errors than naive error analysis. Ta- 
ble [TV] shows the resulting error in the coefficients. Table [V] shows the corresponding error 
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ftdl 




B fuii = B B 


q*a 


2.18 


2.01 


2.15 


2.45 


1.22 



TABLE III. Renormalization scales determined by the Lepage-Mackenzie prescription for the 
axial-vector current A„, for the raw lattice operator O l , and for O^ 11 , are similar. Using this 
prescription for a ratio of matrix elements (as for B B ) is unstable, as described in the text; therefore, 
we choose 2.18 as the scale appropriate for B B . 



in B B . The B B parameter is very insensitive to rather large changes in these parameters. 
Variations of 20% in these parameters change the B B parameter by less than the statistical 
bootstrap errors. It is particularly important that the results are not sensitive to the lattice 
spacing because there are a wide range of possible lattice spacings that could have been 
used: a -1 = 1.94 GeV from the string tension [fi~3|l , a -1 = 2.3 GeV from the p mass [18], and 



a 1 = 2.4 GeV from Upsilon spectroscopy 

To compare the results of B B parameters, in the next section we list our results in terms 



of the one-loop and the two-loop renormalization-group-invariant parameter [35|. We also 
scale B B down to 2.0 GeV for the comparison to some other groups (Sec. |VTD which is 
discussed later. 

To compare at one-loop, we scaled B B and calculated Bb using 

BbM = (j^) B B fa) (4-19) 

B B = a s (i2 2 y( p ^~ 2p ^) B B (p 2 ) (4.20) 

where p is defined in Eq. ( f4.8| ) with the relevant anomalous dimensions listed in Table [TI|. 
Since this is a one-loop calculation, we used 

= ft In (4.21) 

Although a one-loop calculation is traditional, one can also calculate a two- loop 
renormalization-group-invariant B B parameter since the required perturbative calculations 
have been done. 

BbM = {^M ) i 1 + ^ L ~ 2p ^ A ) ) BbM (422) 

Bb = a. W^-^) (l - ^ {p{ >L - 2pUj) BbM (4-23) 

Again p is defined in Eq. (|4.8|) and the relevant anomalous dimensions are listed in Table ||. 
Eq. Q4.15D was used to scale B B and calculate B B to second order. 



In making a comparison to other groups, one can use either B B evaluated at some scale 
or B B . There are disadvantages to both. For the former, either a common scale needs to be 
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+U.UUZ 
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in nnQ 

+U.UUO 


in nnn.8 

+U.UUUO 
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+U.UU4 
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-i-D nns 

TU.UUO 


—0 003 


—0 00Q 
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—0 000Q 


—0 009 










-0.0225 






10% 


-rU.UUUO 

-0.0006 


- 
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-rU.UUUO 

-0.0006 


20% 


_i_n nnnQ 
— n nm 5 


- 




- 
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Z B]V = 


-0.202 






10% 


-pU-UUD 

-0.006 






_ 


- 


-rU.UUO 

-0.006 


20% 


+0.008 
-0.012 










+0.008 
-0.012 








Z Bs = 


-0.137 






10% 








+0.003 
-0.003 


+0.002 
-0.003 


+0.003 
-0.004 


20% 








+0.006 
-0.005 


+0.005 
-0.007 


+0.006 
-0.008 



TABLE IV. The absolute changes from our preferred values of the coefficients Zb l , Zb r , 
Zb N i and Zb s as the parameters q*a, a , m£, and Aqq D , are varied by 10% and 20% first 
individually, and then jointly ("All"), from our preferred values. We do not imply and need not 
assume that the input parameters are known to three significant figures (indeed the coefficients are 
quite insensitive to 20% variations in the values of the parameters); rather, we chose central values 
based on references [p. 37]. 



10% 
20% 



■e = 0.152 
+0.007 
-0.009 
+0.013 
-0.017 



* = 0.154 
+0.007 
-0.009 
+0.013 
-0.017 



t = 0.155 
+0.007 
-0.009 
+0.013 
-0.017 



i = 0.156 
+0.007 
-0.009 
+0.013 
-0.017 



; c = 0.157 
+0.007 
-0.009 
+0.013 
-0.017 



TABLE V. The absolute changes in B B , from Eq. ( p.6Q , due to changes in the coefficients Zb l , 

1 (5 s ) 

Zb r , Zb n , and Zb s as the parameters q*a, a , m£, and Aq^ d are varied jointly by 10% and 20% 
from our preferred values. 
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agreed upon or B B must be scaled. For the latter, the dependence of B B on the choice of n/ 
and A is not negligible; B B can vary by as much as 4% to 5% (see Sec. |V|). This dependence 
is also relevant to using one-loop versus two-loop because the difference in the value of 
between Eqs. ( 4.21 ) and (|4.15|) can vary by as much as 10%. The advantage to comparing 
B B at some scale is that the dependence on n/ and A is less significant (~ 1%, see Table |V|). 
Also, given a value for B B (ml), one can quote a value for B B using either 4 or 5 flavors since 
ml is the boundary between nf = 4 and 5 flavors. These give different constant values of 
B B for the different flavor regimes. One should be explicit about which is quoted. 

Even though the numerical results are for the quenched theory, we use %=5 for /i > m£ 
and n/=4 for fi < ml. There is some evidence from studies of the QCD coupling that the 
effects of omitting dynamical fermions can be modeled by using the correct number of flavors 
in the f3 function (see Sloan Hi] for a review). 



V. SYSTEMATIC ERRORS IN THE MATCHING 

The discussion until now has not revealed any large systematic errors in the perturbative 
matching that could explain the difference between our result and UKQCD's. In this section 
we investigate the systematic error caused by combining the perturbative coefficients for the 
two-point and three-point functions in different ways to form the matching coefficient for 
the B B operator. The UKQCD collaboration found a 20% effect when they changed the 



way they organized their perturbative coefficients . 

We consider three different ways of calculating the coefficients Zb x to consider these 
effects. For convenience, we define the following, where X is one of {L, R, N, S, A}: 

Z x = product of (C f x c zf) 

x (l + ^(d l x \n(q*a) 2 + D l x )) (5.1) 

x 



Lin(Zx) = linearization of (C x c Z d 



An 1 L ' A An 

An 



(d l x M q *a) 2 + D l x )\ (5.2) 



We wish to compare three forms of linearization: "fully linearized" Lin(Z^/Z^), "not lin- 
earized" Zl/Z\ and "partially linearized" Lin(Zi)/Lin(Z J 4) 2 . The UKQCD collaboration 
compared their Lin(Zi)/Lin(Zyi) 2 to Z^/Lin(Z^) 2 when they found their 20% effect in B B . 
Since C^ c is very close to 1, Lin(Z^) is approximately equal to Za- Thus comparing their 
preferred Lin(Z^)/Lin(Z^) 2 to their alternative Z^/Lin(Z^) 2 is essentially the same as com- 
paring Lin(Z i )/Lin(Z^) 2 (partially linearized) to Z L /Z\ (not linearized). 
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To allow a direct comparison with others, our not-linearized results have changed some- 
what from those reported in the conference proceedings [14| and the original preprint of this 
article which calculated the not-linearized result for Zr and as (a^(m£) / a c s (g*)) p °' L Zx 
rather than C{ C Z^. 

In Table [VT] we show the coefficients of the individual B B parameters for the three 
different linearizations described, both with and without tadpole improvement. Table VII 
shows the corresponding value for Bb at both one-loop and at two-loops. The variation 
among the three different linearizations of the non-tadpole-improved coefficients is much 
larger than for the tadpole- improved coefficients. Because there are equal numbers of quarks 
in the numerator and denominator, the individual B B parameters should be independent of 
the wave-function normalization of both the heavy and the light quarks. This implies that 
the coefficients should be independent of tadpole improvement. Tables [VT| and |V11| show 
that this is only true for the fully-linearized quantity, Lin(Z/,/Z^). 

From Table |VJ the overall change in 5^(m^) for the three different linearizations, when 
calculated with the tadpole-improved coefficients, is 20%. However, when calculated from 
non-tadpole-improved perturbative coefficients, Bb{tti\) can change by a much larger fac- 
tor. This suggests that the order-a 2 effects may be large. While these can be treated in 
a variety of ways, we think that they can be treated well or treated poorly. For example, 
the use of tadpole improvement stabilizes the central values and reduces statistical errors. 
The UKQCD collaboration did not use tadpole improvement, which suggests that their per- 
turbative coefficients may be unnecessarily sensitive to their choice of linearization. (Their 
preferred choice is what we call "partially linearized" ; they also considered what we call "not 
linearized".) Their decision not to use tadpole-improvement was forced upon them by the 
way they implemented the light-quark field rotations which were required to remove 0(a) 
corrections to matrix elements [|T5|]. 

We rank the various organizations of perturbation theory in decreasing order of pref- 
erence: fully linearized, not linearized, partially linearized. We discuss, in turn, several 
(related) disadvantages with partial linearizing: larger relative statistical errors, increased 
sensitivity to the value of the lattice coupling constant (via choice of prescription), and 
non-optimal handling of order-a 2 terms. Firstly, due to the larger off-diagonal coefficients 
in the terms of the sum in Eq. ([4.13b ), the numerical result for B B (m\) using non-tadpole- 
improved partially-linearized coefficients (the last row of Table [VT] ) has a larger relative 
statistical error than do the results from the other choices of linearization. 

Secondly, we studied the stability of the results from three groups: the (3 = 6.2 clover- 
static UKQCD simulation ||21|| , the (3 = 6.0 clover-static Gimenez & Martinelli simula- 



tion P5| , and our (3 = 6.0 Wilson-static simulation (both tadpole-improved and not-tadpole 
improved). All three of these groups that have done static Bb simulations used slightly 
different ways of evaluating the perturbative coefficients. We have analyzed all simulation 
data consistently to facilitate comparisons of the results. We compared the linearizations 
for two lattice couplings (a and a v (q*a = 2.18)) and for summing the logarithms (a la 
renormalizat ion-group (RG) techniques) versus not summing the logarithms. These are 
discussed further in Sec. [VI]. We see the same trends in each group's data. Each group's 
partially-linearized result is less stable under variations of a than is their not linearized or 
fully linearized. Their fully-linearized result is close to their not-linearized result; these are 
20% higher than their partial-linearized result. 
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Method 



Zb l Zb r Zbj, 

With Tadpole Improvement 



Zb s 



; n7n +u.uuy 
i.u/u_ 009 

±.UDD_ 0016 
j+0.014 
0.021 



-n n22 +u - UUi 

I +0.002 
-0.004 
h0.003 



-0.008 
-0.012 
-0.021 



-0.008 
-0.013 
-0.025 
-0.015 
-0.021 



B B (ml 



TIM 
0.04 
+0.04 
0.04 
+0.04 
0.04 



Lm(Z x /Z 2 A ) 
(Zx/Zl) 

Lm(Z x )/Lm(Z A ) 2 



-0.202; 



-o.o3i: 



i.oo3: 



-o.o4i: 



0.005 



-0.275_ 037 

n 071 +0.026 
~ u -°' 1 -0.050 



-0.137; 
-0.246; 
-0.251; 



0.98 
0.96 
0.80 



Without Tadpole Improvement 



Lm(Z x /Z 2 A ) 


1 D70 +u - uuy 
i.u/u_ 0i008 


-0 022 +u uul 


-0 202 +u - uus 
u.zuz_ 012 


-0.137j 


-0.006 
-0.007 


0.981-- 


(Zx/Zl) 


1.030lffi 


-0 43+ 004 


PI QQ/l+0-036 
-U.d84_ 0065 


-0.343j 


-0.018 
-0.042 


o.87+g:8S 


Lm(Z x )/Lm(Z A ) 2 


0-8Q2±g5g 


-o n59 +0 005 


U.dZ»_ 092 


-0.358j 


-0.025 
-0.044 





TABLE VI. The effects of different linearizations on the coefficients: The errors on the coef- 

i (5) 

ficients are the statistical errors of varying the parameters q*a, a , mi, and A^ D by 20% from 
our preferred values. The error bars on _Bg(m£) are the bootstrap errors. S^(m^) is the chiral 
extrapolation of the "combine-then-fit" values from Eq. (2.7). 



Method 5(4.33) n f 


one-loop 
A B B 


two-loop 
A B B 


With Tadpole Improvement 


Lm(Z x /Z 2 A ) 0.98(4) 
{Zx/Z 2 A ) 0.96(4) \ 
Un(Z x )/Lm(Z A ) 2 0.80(4) ;> 


175 1.40(6) 
226 1.36(6) 
175 1.37(6) 
226 1.33(6) 
175 1.14(6) 
226 1.11(6) 


175 1.50(6) 
246 1.46(6) 
175 1.47(6) 
246 1.43(6) 
175 1.23(6) 
246 1.19(6) 


Without Tadpole Improvement 


Lm(Z x /Z 2 A ) 0.98(4) |j 
(Zx/Zl) 0.87(4) \ 
Lm(Z x )/Un{Z A ) 2 0.49(4) | 


175 1.40(6) 
226 1.36(6) 
175 1.24(6) 
226 1.21(6) 
175 0.70(6) 
226 0.68(6) 


175 1.50(6) 
246 1.46(6) 
175 1.34(6) 
246 1.30(6) 
175 0.75(6) 
246 0.73(6) 



TABLE VII. From the B 



B 



ml) result extracted by Monte Carlo, listed in Table VI, we calcu- 
lated a Bb with both 4 and 5 flavors (see text). The Lin(Z x /Z A ) results are our preferred values. 
As mentioned in the text, Bb varies with nf and A^ n f^ as well as with loop-order. 
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Thirdly, we believe that partial linearization does a poor job of organizing higher-order 
terms. The treatment of 0(a 2 ) terms in partially-linearized coefficients causes the low values 
seen by all groups by linearizing some terms but not the whole ratio. We prefer the fully- 
linearized method because it removes all of these 0(a 2 ) terms (as in a Taylor expansion) by 
linearizing the whole ratio. Fully- or not-linearizing the coefficients treats the 0(a 2 ) terms 
more appropriately than does partially linearizing. 

Our preferred choice of linearization (full) can also be motivated by the non-perturbative 
renormalization method, introduced by the Rome- Southampton group [|6||. The non- 



perturbative renormalization method for B B parameter would be very similar to that used 
to obtain the renormalization constants for the kaon B parameter pT7| , fi8|| , in which all the 
factors of the lattice wave function normalization of the quarks cancel explicitly for the B 
parameter. In perturbation theory, this corresponds to our preferred full linearization. The 
non-perturbative method only determines the lattice part of the renormalization factor; a 
choice of linearization would still have to be made for the continuum factor. However, the 



continuum factor can and should be calculated to next to leading order |46j, making it less 
sensitive to the different choices of linearizations. 

In summary, our preference for the treatment of the coefficients is to fully linearize the 
ratio (in the notation of this section, Z Bl is Lin(Z^/Z^)). This gives a result which has 
no order-a 2 terms, which is insensitive to the inclusion of tadpole improvement and to the 
wave-function normalization model by allowing explicit cancelations, and which reduces the 
statistical errors in B B . The quantitative consequences of our choice are discussed in the 
following section where we compare the results of different groups. 

Just as the numerical value of B B is stable because of cancelations of correlated fluc- 
tuations in numerator and denominator, we have argued that so too are its perturbative 
corrections when fully linearized. The fully- linearized perturbatively-calculated coefficients 
for B B are likely more reliable than those for the product B B f B , the quantity which is 
required in the analysis of B°—B° mixing experiments. In the Appendix, we discuss our 
recommendation for how to linearize the product B B f B . 



VI. WORLD COMPARISON 

In Table |VIII| (|IXD, we show a collection of results from several groups scaled to give 
B B (ml), B B (2.0GeV), and the one-loop (two-loop) renormalization-group-invariant Bb pa- 
rameter. Results from both static and relativistic-quark simulations are shown. The sim- 
ulations using relativistic heavy Wilson quarks PJ^, p3| , p][ ] calculate the B B parameter for 
quark masses around charm and extrapolate up to the physical mass, using a fit model of 
the form 

Bb = B ° B + M (0) 

The value of B° B should be the same as the static theory result. (It is better to do a 
combined analysis of relativistic and static quarks to obtain a value for B B .) We call B B 
the "extrapolated-static" value. 

Tables [VIII] and |X| show that values for B B obtained from Wilson action simulations 



are basically consistent; the small differences can be explained by small lattice-spacing and 
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Method 



^2 

Ref. (GeV) B(u, 2 ) 



A 

(MeV) 



5(2.0) 



one-loop 
5(4.33) 



B B 



Static-Clover pj| 



Static-Clover gl] 



Static-Clover g§ 



Static-Clover (45 ] 



6.2 m b =5.0 0.69(4) 

6.2 m b =5.0 0.81(4) 

6.0 m b =5.0 0.54(4) 

6.0 m b =5.0 0.76(5) 



130 
200 



0.75(4) 0.70(4) 



151 

200 



0.82(5) 0.77(5) 



1.02(6) 
0.98(6) 



5 


130 






1.19(6) 


4 


200 


0.87(4) 


0.82(4) 


1.14(6) 


5 


151 






0.79(6) 


4 


200 


0.59(4) 


0.55(4) 


0.77(6) 



1.11(7) 
1.08(7) 



Static- Wilson 



this 
work 



6.0 m*=4.33 0.98(4) 



175 
226 



1.05(4) 



0.98(4) 



1.40(6) 
1.36(6) 



Extrap. Static g 
Extrap. Static 



5.7-6.3 ^=2.0 1.04(5) 
6.4 u.=3.7 0.90(5) 



200 
226 



1.04(5) 



0.97(5) 
0.97(5) 



200 
200 



0.94(5) 
0.95(5) 



0.89(5) 
0.89(5) 



1.36(7) 
1.34(6) 



1.21(7) 
1.25(7) 



Wilson-Wilson @ 
Wilson-Wilson @| 



5.7-6.3 fi=2.0 0.96(6) 
6.1 ^=2.0 1.01(15) 



Wilson- Wilson @ 6.1 m b =5.0 0.895(47) 



Wilson- Wilson @ 



Wilson- Wilson @ 



6.3 m b =5.0 0.840(60) 

6.4 ^=3.7 0.86(5) 



200 
226 



0.96(6) 



0.90(6) 
0.89(6) 



200 
226 



1.01(15) 



0.94(13) 
0.94(14) 



239 
239 
183 



0.96(5) 
0.98(5) 



0.90(5) 
0.91(5) 



246 
246 
189 



0.90(6) 
0.92(6) 



0.85(6) 
0.85(6) 



200 
200 



0.90(5) 
0.91(5) 



0.85(5) 
0.85(5) 



1.25(8) 
1.24(8) 



1.32(20) 
1.30(19) 



1.21(6) 
1.25(7) 
1.29(7) 



1.14(8) 
1.17(8) 
1.20(9) 



1.16(7) 
1.19(7) 



Sum Rule 



51 



m b =4.6 1.00(15) 



175 
227 



1.08(16) 1.00(15) 



1.43(22) 
1.39(21) 



TABLE VIII. The authors' numbers, quoted at the listed value for /Z2, have been scaled using 
Eq. ( [4.19D to /u=2.0 GeV and to m£=4.33 GeV. The slanted numbers are those that the cited 
authors quote. We calculated B B (ml) in the Static- Wilson case and then scaled it to 2.0 GeV 
using n,f=4 and calculated a Bb with both 4 and 5 flavors (see text). The value quoted by this 
work uses the fully-linearized tadpole-improved coefficients. The JLQCD collaboration cite their 
A's as n/=0 values. When Abada et al. quotes a Bb for the Wilson quarks, they use ^/=0. We 
scaled both groups' results using both n/=0 and nf=4. 
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Method 



^2 

Ref. p (GeV) B(fi 2 ) 



rn 



A 

(MeV) 



5(2.0) 



two-loop 
5(4.33) 



B B 



Static-Clover fl2l 



Static-Clover |21] 



Static-Clover @ 



Static-Clover g5] 



6.2 m fe =5.0 0.69(4J 

6.2 m b =5.0 0.81(4) 

6.0 m fe =5.0 0.54(4j 

6.0 m b =5.0 0.76(5) 



130 
200 



0.74(5) 0.70(4) 



136 
200 



0.77(5) 0.81(5) 



1.09(6) 
1.05(6) 



5 


130 






1.27(6) 


4 


200 


0.86(4) 


0.81(4) 


1.23(6) 


5 


136 






0.86(6) 


4 


200 


0.58(4) 


0.54(4) 


0.82(6) 



1.21(8) 
1.16(8) 



Static- Wilson 



this 
work 



6.0 m£=4.33 0.98(4) 



175 
246 



1.04(4) 



0.98(4) 



1.50(6) 
1.46(6) 



Extrap. Static g 
Extrap. Static 



5.7-6.3 /z=2.0 1.04(5) 
6.4 /x=3.7 0.90(5) 



200 
246 



1.04(5) 



0.98(5) 
0.98(5) 



200 
200 



0.93(5) 
0.94(5) 



0.89(5) 
0.89(5) 



1.49(7) 
1.46(7) 



1.29(7) 
1.35(7) 



Wilson-Wilson g 
Wilson-Wilson f|| 



5.7-6.3 /z=2.0 0.96(6) 
6.1 /z=2.0 1.01(15) 



Wilson- Wilson @ 6.1 m fe =5.0 0.895(47) 



Wilson- Wilson @ 



Wilson- Wilson 



6.3 m 6 =5.0 0.840(60) 

6.4 ^=3.7 O.S6(5j 



200 
246 



0.96(6) 



0.91(6) 
0.90(6) 



200 
246 



1.01(15) 



0.96(14) 
0.95(14) 



239 
239 
183 



0.94(5) 
0.96(5) 



0.90(5) 
0.90(5) 



246 
246 
189 



0.88(6) 
0.90(6) 



0.85(6) 
0.85(6) 



200 
200 



0.89(5) 
0.90(5) 



0.85(5) 
0.85(5) 



1.37(9) 
1.35(9) 



1.44(21) 
1.42(21) 



1.29(7) 
1.35(7) 
1.36(7) 



1.21(9) 
1.26(9) 
1.30(9) 



1.24(7) 
1.29(7) 



Sum Rule 



m 



m b =4.6 1.00(15) 



175 
227 



1.07(16) 1.00(15) 



1.54(23) 
1.50(22) 



TABLE IX. This table repeats the analysis in Table VIII, using the two-loop renormalization 
group invariant B B parameter. 
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5^(4.33 GeV) 


UKQCD gj| G&M H TAD NO-TAD 


s 5 

fl (M3) 


0.84(5) 0.85(4) 0.97(4) 0.97(4) 
0.83(5) 0.85(3) 0.97(4) 0.97(4) 
0.75(5) 0.84(4) 0.96(4) 0.96(4) 
0.77(5) 0.84(3) 0.96(4) 0.96(4) 


x a 

Oty 


E 5 

111 (Ml) y — 


0.85(5) 0.84(3) 0.95(4) 0.81(4) 
0.82(5) 0.82(3) 0.96(4) 0.87(4) 
U.7o(oJ 0.83(3J 0.94(4J 0.80(4J 
0.76(5) 0.81(3) 0.95(4) 0.86(4) 


e a 


E 5 
i a y 


0.72(5) 0.70(3) 0.75(4) 0.30(3) 
0.62(4) 0.62(3) 0.80(4) 0.49(4) 
0.60(4) 0.68(3) 0.73(4) 0.27(3) 
0.54(4) 0.61(3) 0.78(4) 0.47(4) 


pi (M2) ?— 

E a 

OLy 



TABLE X. Comparison between the fit-then-combine (Eq. ([2.6;)) analysis for £?(m&) of the 
three groups' data. These numbers are for m& = 4.33 GeV, q* = 2.18a _1 , and nj = 5. "fl" is 
fully-linearized, "nl" is not-linearized, and "pi" is partially-linearized. "Ml," "M2," and "M3" refer 
to the notation of Gimenez & Martinelli and Wittig J5^]. We list both our tadpole-improved 
(tad) and our non-tadpole-improved (no-tad) results. The errors are roughly estimated from 
statistical errors on the raw Bx values and approximate errors on the coefficients. See the text for 
comments on a and et v . 



finite-volume effects. Our result is consistent with that of Bernard and Soni, as reported by 
Soni || , for the extrapolated static Wilson fermions. 

Since the original appearance of the preprint for this article ([hep-lat/96 100261 version 
1), data has been made available which allows a more detailed comparison between ourselves 
(on the high end of the world results) and others (on the low end). Firstly, we have added 
the updated numbers from Gimenez & Martinelli H5[] to Tables |V111| and [TX| Secondly, we 
note that Wittig [0] has a nice review on the subject of leptonic decays of lattice heavy 
quarks, in which he compares the results of UKQCD pi |, Gimenez & Martinelli and 
the preprint of this article. 

In his Sec. 4.2, Wittig offers Table 9 for comparison, using our non-tadpole-improved 
results. We find that the tadpole-improved Wilson results improve the non-tadpole-improved 
results, so we prefer to compare their clover-improved results to our tadpole- improved results. 
Our analogous comparison results in the numbers listed in Table [X|. 

In the comparison, we use rif — 5 and Aqq D = 0.175 GeV which result in a c s ont (rrib = 
4.33 GeV) = 0.21. We also use our two-loop A v a = 0.169 to scale a l $ tt (q* = 2. 18a -1 ) = 0.18. 
Both a s are with two-loops from Eq. (fj.l5Q . We also compare using a = 6 / (4ti (3uq) , which 
is 0.132 for the UKQCD collaboration 0.1458 for Gimenez & Martinelli |5J, and 0.198 
for us. (For each group, we used u = l/8K crit to calculate a.) In addition, since the 
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original UKQCD results! do not sum the logarithms (a la RG techniques), Table |X| lists 
both summing logs (E) and not summing the logs (S). 

Rather than calculate a q*a (Eq. ( f4.17| )) and a A v a (Eq. ( [4.14| )) for the clover action, 
we used our values. Since a v is a function of (q*a/ A v a), a v (q*) is the same for all three 
groups. We note that a -1 = 2.9 was used for UKQCD and a -1 = 2.1 was used for both 
Gimenez & Martinelli and ourselves. Since q*a was chosen to be 2.18 for all three groups, 
the scales in the comparison of Table [X] are different. This is the reason that the UKQCD 
2 results differ from their £ results. The $ results are more sensitive to the scale of the 
perturbative matching. 

Though not listed in the table, we are able to reproduce the results of both UKQCD f2T| 
and Gimenez & Martinelli [45] for m;, = 5.0GeV, /i = a~ 1 , nj = 4, and Aqq D = 0.200 GeV 
when we tailor the respective calculations according to the method presented in each paperi. 
Also, we agree with the results of Wittig [52] for our Tt-a v entries when we use his param- 
eters. 

Both UKQCD's and Gimenez's & Martinelli's quoted values for the static B B are lower 
than all of the other results. One possible reason for these low results is that they used 
the clover action for the light quarks, which does not have corrections to the continuum 
limit that are linear in the lattice spacing, whereas the standard Wilson fermion action does 
have such artifact terms. However, the Wilson results are stable over four different lattice 
spacings, which implies that the lattice artifact terms alone cannot account for the difference 
between the clover results and the Wilson numbers. 

Table [X] shows that the not-linearized (and fully-linearized) static clover results for B B 
are larger than the partially-linearized results, as is discussed in the original papers. The 
clover-static results that use the not-linearized matching are in better agreement, though 
still low, with the results from simulations which use relativistic heavy quarks to simulate the 
b quark (see Table |V111| ). All the published data flP5|] on calculating B B using relativistic 
heavy quarks favor a negative value of B B in Eq. ( |6.1|) . For consistency, the static value of 
B B should be higher than the value of B B extrapolated to the b quark mass. This is true 
for our result and favors the higher clover-static results. 

The various choices made in the calculation have non-negligible effects. One can choose 
which action to use (Wilson vs clover), whether or not to tadpole-improve, and which 
linearization method to use. The choice between our tadpole-improved Wilson-static action 
and the non-tadpole-improved clover-static action has a 15% effect in both the fully- and not- 
linearized (S-a^) cases. This is a 20% effect for the partially-linearized case. In addition, 
tadpole- improvement stabilizes the Wilson-static results to the extent that one can make a 



3 UKQCD did investigate the use of renormalization group improved perturbation theory, but they 
did not use it to calculate their final results. 



To reproduce UKQCD's |21| 0.69(4) and 0.81 (the latter is our conversion of their quoted B = 
1.19), do not sum the logs, use /za=l, and do not include the cross term, (U T ) LS , in the coefficient 



of Ol- To reproduce Gimenez' & Martinelli's [45] Table 3, sum the logs and include the cross term, 
but use fia=l, even for the a v (q*a = 2.\8) case. 
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better comparison of different linearizations between tadpole- Wilson-static and clover-static 
than between non-tadpole- Wilson-static and clover-static. Finally, there is a 20% effect due 
to choice of linearization for either action. This linearization effect is at least as large as the 
effect due to choice of action. For reasons given in Sec. [V|, our favorite choice of linearization 
is the fully-linearized treatment. 

A similar trend can be seen in each group's results: partially-linearized values are smaller 
and less stable than either not-linearized or fully-linearized values. This is due to 0(a 2 ) terms 
which may or may not cancel to varying degrees. The partially-linearized treatment only 
linearizes part of the ratio which causes its value to be misleadingly low. The not-linearized 
and fully-linearized treatments are better because they do not do this. The fully-linearized 
treatment is preferred because it treats 0(a 2 ) terms uniformly by removing them (as one 
does in an expansion). 



VII. CONCLUSION 

Our primary result from this tadpole-improved j3 = 6.0 Wilson-static calculation is 
B B {ml) = 0.98i|i? 8 , where the errors are statistical (bootstrap) and systematic, respec- 
tively. The overall systematic error is obtained in quadrature from the following: 1J from 
the choice of fit-range, from the parameter-dependence of the perturbative-calculated 
mixing coefficients, and i?g due to the the choice of linearization of the coefficients, as 
was discussed in Sec. |V|. The unusual asymmetry of the latter systematic error reflects 
our preference for a particular choice of linearization ("full"). Our second favorite choice 
("not-linearized") results in a central value of 0.96. We quote a very conservative systematic 
error to encompass our least favorite choice (0.80 from "partial linearization") even though 
we have argued against this choice. Systematic errors from finite lattice spacing and from 
quenching are not estimated. 

Tables [VII]] and fX] show that values for B B obtained from Wilson action simulations 



are basically consistent; the small differences can be explained by small lattice-spacing and 
finite- volume effects. The simulations all favor a negative value of B l B in Eq. (|6.I|) ||. For 
consistency, this implies that the static value of B B should be higher than the value of B B 
extrapolated to the b quark mass. Our number is on the high end of the comparison in 
Table |V111| and is consistent with that of Bernard and Soni who use extrapolated static 
Wilson fermions. 

In Sec. [V] we investigated the effect of changing the way the four-fermion operator renor- 
malization and the axial-current renormalization were combined to form the matching co- 
efficient for the B B parameter. We presented arguments that suggested that our preferred 
way of organizing the continuum-to-lattice matching (full linearization) was superior to any 
other method we considered. We also showed that making a different choice could lower 
the result by as much as 20%. Besides the linearizations, Table [X| shows a 15% difference 
due to choice of action between our tadpole- improved (3 = 6.0 Wilson-static and the non- 
tadpole- improved (3 = 6.0 and 6.2 clover-static results in the fully- and not-linearized cases. 
(The Wilson results are at the high end of the world data and the clover results are at the 
low end.) Partial-linearization leaves a 20% effect due to choice of action. The effect due to 
choice of linearization is at least as large as the effect due to choice of action. 
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Although all organizations of perturbation theory at one-loop are theoretically equal, 
some are more equal than others! Fully linearizing gives a result which has no order-a 2 terms 
and which is insensitive to the inclusion of tadpole improvement and to the wave-function 
normalization model by allowing explicit cancelations (which reduces the statistical errors 
in B B ). 

In our perturbative-matching procedure we included next-to-leading order log terms 
and organized the perturbative matching in a way that we believe reduces higher-order 
corrections. Also we used the automatic scale-setting procedure of Lepage and Mackenzie 
to find the "best" scale to use in the lattice-to-continuum matching. The agreement of 
our results with relativistic heavy quark results supports our procedure. Our conclusion 
is that for the Wilson-static case, the use of tadpole improvement and of a fully-linearized 
treatment of the mixing coefficients is preferred. Of course, this may become less important 
numerically with increased coupling and/or improved actions; however, we still recommend 
the procedure. 

Although sensible things can be done to reduce the effects of higher-order perturbative 
corrections in the lattice-to-continuum matching, this will remain the dominant uncertainty 
in the calculation of B B in the static theory. In principle, "all" that is required is a calculation 
of the two-loop anomalous dimension of the Ol and operators on the lattice. Although 
this calculation is very difficult, new developments in lattice perturbation theory for Wilson 
quarks |53| and a new stochastic way of doing lattice perturbation theory |54[] may make 
these calculations more tractable in the future. A more immediate solution would be to use 



the numerical renormalization technique, developed by the Rome- Southampton group [46 



which has already been used to determine the lattice perturbative coefficients for static 



f B |[48|| , for the kaon parameter |47]], and for other important quantities. 

The relative consistency of the Wilson B B results motivates a large study using both 
relativistic and static quarks in the same simulations to constrain the interpolation to the 
B mass. To constrain the systematic errors, the results of simulations with different lattice 
spacings and volumes should be combined to take the continuum limit. This kind of study 
will also help to control the perturbative-matching errors, as the effects of the higher-order 
perturbative terms are reduced as the continuum limit is taken. (A nice example of this 
for the effects of different renormalization prescriptions on light-quark decay constants has 



been given by the GF11 group |5"5| .) 

Once mixing in the B®—B° s system has been measured experimentally, the results can 
be combined with data from B°—B° mixing experiments to calculate the V ts /V t d ratio of 
CKM matrix elements. The advantage of calculating this ratio is that various uncertain 
standard-model factors cancel. However, a value of BbJbj ^b/b is required. As there are 
a large number of lattice results on the calculation of fs 3 / fB H31 > nere we concentrate on 
the ratio BbJBb- 

Using a fit model which is linear in the quark mass, we obtain B B JB B = 0.99+1(1). (The 
first error is statistical (bootstrap) and the second is the standard deviation of the fitted 
value for "reasonable choices" of fit range.) Even though the ratio B B JB B is determined 
quite precisely, it is not resolved whether B Bs is greater than or less than B B since the 
B B parameter is found to depend weakly on the quark mass. Other groups pl| , |50| , |56| have 
reported similar findings. Most lattice simulations have found that fs s is between ten 
and twenty percent larger than f B [TI|. Bernard et al. ]5E[ have extracted the ratio of 
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Bsf^jBsfs directly by doing individual fits to the three-point function in relativistic quark 
simulations. This is a promising approach for relativistic heavy quarks. We did not try it 
because of concerns about the signal-to-noise ratio and about the size of the perturbative 
coefficients in the static theory. 

Our result also contains an unknown systematic error due to quenching. Quenched 
chiral perturbation theory predicts the effects of quenching to be small for B B |57| , |58| ; this 
conclusion was confirmed by Bernard and Soni |2j who calculated B B in both quenched and 
dynamical simulations. In Soni's review || of the lattice calculation of weak matrix elements 
at the Lattice '95 conference, he quotes a value of B B {2 GeV) = 1.0 ±0.15 (90% confidence) 
as his best estimate of the B B parameter. Our result, 5 B (2GeV) = 1.05i 4 jl?9, * s consistent 
with this value and with the vacuum-saturation-approximation value, 1. 



This work is supported in part by the U.S. Department of Energy under grant numbers 
DE-FG05-84ER40154 and DE-FC02-91ER75661, and by the University of Kentucky Center 
for Computational Sciences. The computations were carried out at NERSC. 



In the analysis of B°—B° mixing experiments the value of B B f B is required. Here we 
discuss the linearization options for combining B B and fg from a variety of linearizations 
of these quantities. If a not-linearized B B is multiplied with a not-linearized (/b) 2 , then 
the only order-a 2 effects which remain are due specifically to not linearizing the numerator 
of B B . We estimate this effect to be on the order of 10%. If one multiplies a partially- 
linearized B B with a linearized Lin(ZA) 2 , then there should be no order-a 2 effects due to 
the product. However, if one mixes a linearized with a not-linearized B B and then there 
can be terms of almost 20%. Although the difference between Za and Lin(Z^) is smaller 
than 5%, the difference between Lin(Z^) 2 and Lin(Z^) is just over 15%. The practical 
drawback of using a B B which is not linearized or is partially linearized is that there are 
order-a 2 terms present which may or may not cancel when the B B is combined with an 

The practical drawback to using the fully-linearized B B is linearizing the product B B f B . 
This is easily remedied. The fully- linearized B B , Bfi, essentially has the form 



where the Br, Bn, and B$ can be included by adjusting the values of A and C appropriately. 
When this is combined with the square of the linearized /, 
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APPENDIX: LINEARIZATION STRATEGY FOR B B F% 




(Al) 




(A2) 



it would be convenient to get a linearized result with no order-a 2 terms: 




(A3) 
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Since 



[1 + aA) [l + a 



B 



[1 + aA + aB) 



1 + aA 

this is straightforward to accomplish. The product of Bfi with the linearized square of 

t \ 



(A4) 



fL 



lin 



l + a c 



D 



\ 



+ a 1 



E 



B 



r 



(AS) 



r L aw value can 



gives Eq. ( |A3| ) with no order-a 2 terms due to coefficient multiplication. Our B 
be read from the first row of Table |. 

While the product of the partially-linearized B B with the linearized Jb also does not have 
any order-a 2 terms due to coefficient multiplication, the partially-linearized B B by itself has 
order-a 2 terms which are on the order of 18% (See Tables |VI] and [VI I|) . The advantage of 
our method is that all three quantities BB(ml), /e(m^) and 5s/^(m^) have no order-a 2 
terms due to coefficient multiplication, and that Ss(m^) is stable against the inclusion of 
tadpole improvement and the choice of wave-function normalization. 
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FIG. 1. Effective mass m(t + 1/2) = lnC(t)/C(t + 1) from the LS (local sink, smeared source) 
two-point correlation function C(t). The source was smeared with an optimal smearing function 
produced by the most [12| algorithm which was designed to eliminate excited-state contamination. 
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FIG. 2. Same as for Fig. |l] but for the SL (smeared sink, local source) two-point correlation 
function. The same optimal smearing function is used to eliminate excited state contamination, 
but statistical errors are larger since the source is (necessarily) a delta function. 
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FIG. 3. Schematic diagram of the quark flow for the three-point correlation function of 
Eq. ( |2.1| ), The "targets" are intended to represent the smearing of the light quark relative to 
the static quark. The static quarks are restricted to the spatial origin. 
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FIG. 4. Raw £> parameter for the Ol operator from Eq. (|2.5|). 
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FIG. 5. The dependence on the fitted raw Bl parameter on the choice of ii, the (fixed) time 
position of one interpolating field, and on the fit range t 2 mm — 1 2 max of the other. Clustered points 
have different £2, max- All fits take into account correlations in t 2 , and are not displayed if the naive 
quality of fit Q does not exceed 0.2. 
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FIG. 6. Same as for Fig. ^| but for the On operator. 
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FIG. 7. Same as for Fig. || but for the On operator. 
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FIG. 8. Same as for Fig. ||] but for the Os operator. Note the normalization as explained in 
Table I 
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FIG. 9. The ratio of (the linear combination of) three-point functions to two-point functions 
which approaches B B for large Euclidean times. 
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FIG. 10. Same as for Fig. [| but for the B B parameter itself. 
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